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In a recent Comment [T], Bizon and Rostworowski 
present a few criticisms of our recent Letter [2]. In par¬ 
ticular, they present three arguments: (1) that their own 
evolutions of the two-mode initial data we had studied 
collapse to a black hole around a time t ss 1080 whereas 
we found no collapse for times up to roughly t = 1500, 
(2) that our two-timescale framework (TTF) u cannot he 
even used to infer stability ,” and (3) that our 11 claims 
about stability islands” may not be correct. 

We have studied evolutions of this initial data with res¬ 
olutions higher than we had originally. As displayed in 
Fig. □ we do find that higher resolutions display higher 
concentrations of energy, but we nevertheless have still 
not observed collapse to a black hole. We note that our 
code demonstrates convergence to a unique solution even 
at the late times in question (i.e., around t ~ 1080). We 
also emphasize that we have demonstrated that the solu¬ 
tion to which our numerical results converge is in fact a 
solution of the scalar anti-de Sitter (AdS) system we seek 
to solve. This latter property ( consistency) is demon¬ 
strated with tests that the constraint residuals and mass 
loss converge to zero (see our Supplementary Material). 
And so, while it is possible that the continuum solution 
does in fact collapse around t « 1080 (and indeed, the 
work of [3j does see collapse in a code for which both con¬ 
vergence and consistency are verified), we cannot confirm 
this (in time for the publication of [I]) and hence their 
claim (1). 

However even if collapse does take place at a time t ss 
1080, the main claims of our paper still stand. We could 
not know if collapse occurred after the time which we ran 
our code, and so collapse for some time soon after t = 
1500 was always a possibility. And so if instead collapse 
occurs at f ~ 1080, there is no change to our claims. 
Within that time one still finds both direct and inverse 
cascades which is now confirmed by the Comment. 

Another important point is that, were we to decrease 
the initial amplitude e, we would certainly observe any 
possible collapse pushed to times later than t ss 1080. 

Regarding their claims (2) and (3), there seems to be 
some misunderstanding of what we tried to communicate 
in our Letter. Determining whether some scalar pertur¬ 
bation of AdS with arbitrarily small amplitude e collapses 
to a black hole is not possible with numerical evolutions 
which require some finite e, finite resolution, and finite 
evolution time t. 

More particularly regarding their claim (2), we agree 


that the TTF cannot be used to infer stability beyond 
times oc e -2 . Indeed, we stressed that we have carried 
our TTF analysis to 0(e 3 ) so that one can only trust 
predictions within times scaling as e~ 2 . However, for 
shorter times, we disagree with the claims of the Com¬ 
ment that our truncation to j max = 47 “does not suffice 
to capture the dynamics of the turbulent cascade” for 2- 
mode initial data. As is clear from Fig. 4 of [2], the vast 
majority of the energy remains in the lowest modes of 
the system during evolution, and the dynamics of these 
modes (in particular, the recurrence times) are very well- 
captured by the TTF evolution with j max = 47. We have 
clearly stated the limitations of our mode-truncation in 
the Supplementary Material; in particular, we explained 
that this leads to the discrepancy with the full numeri¬ 
cal GR curve in Fig. 3 of [23 because the high-j modes 
are highly peaked about the origin (despite carrying very 
small amounts of energy). We note that the TTF has 
been verified in Ref. [3] and even the authors of the Com¬ 
ment have now used the method [5j. 

With regard to claim (3), we did not use the phrase 
“two-mode stability islands” in our Letter. In later 
work |6] we referred to “stability islands,” but in that 
case we were referring to quasi-periodic solutions. 

Their Comment provides independent confirmation of 
our main claims: (i) the presence of both direct and in¬ 
verse energy cascades (a fact that was missed in the orig¬ 
inal perturbative analysis), and (ii) the validity of the 
TTF. 
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FIG. 1. The behavior of the scalar field at the origin for the 
evolution of equal-energy, two-mode initial data for e = 0.09 
with increasing resolution. 









